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Abstract 

CN ! We calculate cross sections for gg — >• QQQQ in the high-energy approximation in the mixed 

j^— ^ I (longitudinal momentum fraction, impact parameter) and momentum space representations. 

CN ' Besides the total cross section as a function of subsystem energy also differential distributions 

H ' (in quark rapidity, transverse momentum, QQ, QQ invariant mass) are presented. The clemen- 

cy ■ 

tary cross section is used to calculate production of {cc){cc) in single-parton scattering (SPS) 

in proton-proton collisions. We present integrated cross section as a function of proton-proton 

center of mass energy as well as differential distribution in Mr^ii-jr^gj. The results are compared 

with corresponding results for double-parton scattering (DPS) discussed recently in the liter- 

' (— i '' ature. We find that the considered SPS contribution to {cc){cc) production is at high energy 

0-c {^/s > 5 TeV) much smaller than that for DPS contribution. 
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I. INTRODUCTION 



The cross section for cc production in proton-proton or proton-antiproton collisions 
at high energy is quite large [l|, |2| because the gluon-gluon luminosity grows quickly with 
energy. We have shown recently that the cross section for production of two pairs of cc 
in double-parton scattering (DPS) grows even faster with incident center-of-mass energy 
[3| and becomes very large at large energies. In order to verify the DPS contribution 
a single-parton scattering (SPS) contribution has to be evaluated. This was not done 
so far in the literature as it requires calculation of 2 — )■ 4 subprocesses. At LHC a 
measurement of the two-pairs of cc production should be possible. This could be a good 
test of methods of calculating higher-order QCD corrections. 

It is the aim of the present paper to calculate contribution of single-parton scattering 
to the inclusive pp — )■ {cc){cc)X cross section. In the present paper we shall use high- 
energy approximation in calculating elementary gg — )■ {cc){cc) cross section. At low 
incident energy and/or low cccc invariant mass production a carefull treatment of the 
threshold effects is required. The elementary 2 — )■ 4 cross section is convoluted next 
with gluon distribution functions. The result is compared with that for DPS presented 
recently |3|. A prospects how to disentangle SPS and DPS contributions will be discussed 
in the Result section. 



II. THEORETICAL FRAMEWORK 
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FIG. 1: One-gluon exchange amplitude for the gg — )■ {cc){cc) process. 



A. gg -^ QQQQ 



In the high energy limit of the gluon-gluon subprocess, the gg — )■ QQQQ amplitude 
will be dominated by the t-channel gluon exchange. Furthermore, if the QQ-pairs are 
produced in the respective fragmentation regions of the incoming gluons an intuitively 
appealing approach based on light-cone perturbation theory is possible. Namely we are 
looking for the cross section of excitation of the QQ-Fock states of the colliding gluons. 



We base our calculations on previous works [4, |5[, where the total cross section for 
the process gN — )■ QQX and gA — )■ QQX has been obtained, as well as on |6|, where 
more detailed kinematical distributions can be found. 

Let us briefly recapitulate the derivation of the key formulas. We start by writing out 
the lightcone Fock-state expansion of the incoming, physical, gluon, takiiig into account 
only the heavy-quark-antiquark fluctuation in the two-body sector (see, |5|): 

\g^{b)) = ,/T^^\g^,J + I d'rdz^{z,r)\[QQr,;z,r) . (1) 

Here, quark and antiquark in the gluon carry fractions z,l — z oi the gluon's large light- 
cone plus-momentum and are seperated by a distance r in the impact parameter plane. 
Notice that the quark-antiquark pair before the interaction is not a color dipole, but 
carries a color-octet charge. Infrared safety of the relevant cross section is a consequence 
of the fact, that arbitrarily long-wavelength gluons cannot induce the g — )■ QQ transition 
and decouple. 

The normalized color-states of the quark-antiquark system in the color-octet and 
color-singlet states are, respectively: 

iiQQrs) = V2{n] UQ') , i[gg]i) = ^^} UQ') ■ (2) 

VNc 

A similar Fock-space expansion as ([I]) can be written for the second gluon participating 
in the interaction, which will have a large momentum in the light-cone minus-direction. 
The interaction between the right- and left moving parton systems is then mediated 
by the gluon exchange which acts like a helicity conserving potential [TJ between partons 
i and j. 

V{b + h - s,) = (-.)^ j p^ eM^ib + h- s,)q] T^^T^. (3) 

Here b is the impact parameter between the two initial gluons, 6j and Sj are the 
impact parameter distances of Q and Q relative to their respective parent gluon which 
are conserved during the interaction. The matrices T^ are the generators of color SU{Nc) 
acting on parton i in the relevant representation. Notice, that the gluon mass parameter 
fiG is not needed to make the cross section finite, as mentioned above there are no 
singularities associated with the small-q behaviour of the gluon propagator. As a physical 
parameter it enforces a finite propagation radius of gluons in the transverse plane, as 
is in fact enforced by confinement. In practice, for the problem at hand its precise 
value is unimportant: as long as /i|. ^ uiq our results are practically independent of 
fiG- The relevant Feynman diagrams for our process are shown in Figs.([T]) and ([2]), and 
the corresponding scattering amplitude for gg — )■ {QQ){QQ) in the high-energy limit of 
interest takes the form: 



[S{b) - 1] \g^ ® /; b) = -^^ [ f^^ \ J^ e"^^^^ T^^] ® \ ^ e^^^'^ TJ:\g^) 



(4) 



Let us now turn to one of the factors in the square brackets (the so-called "impact- 
factors"). The quark in the gluon is located at the distance bq = {1 — z)r, and the 
antiquark at bq = —zr, where r = bq — bq. In fact we are interested in the respective 
two-body Fock-state components orthogonal to the physical gluon, and the relevant piece 
of the amplitude for the excitation g — )■ QQ is then given by (see also (Sj): 



excitation 



k 

= I dzd'r^{z,r){[e'^'-^^'''l ~e-"''1]-^6a,\[QQ],;z,r) 

+ [e^(i-)r-g + e-^^^-2]lz/,,,|[QQ]-z,r)} . (5) 

Here b is the color index of the t-channel gluon which carries the transverse momentum 
q. Clearly, the impact factor vanishes for q — t- 0. A completely analogous expression can 
be written for the other factor in square brackets, and the full amplitude is obtained as: 

A(^ V ^ QQQQ; &) = -^— [ -^-^^^\q, g^ ^\-q, /) • (6) 

The two impact factors correspond to the upper and lower gluons, respectively. The total 
cross section is then obtained after integrating the squared amplitude over the impact 
parameter and averaging over initial gluon colors 



crtot = ^^, \ ^^^ ^y'rf^b|A((?v -> QQQQ; b)\' = 4alj 



[q' + f^l? 



(7) 



Similar impact factor representations for related QED problems are known for a long 
time J9|. Introducing the short-hand notation 

F{qr) = 2 — exp(iqrr) — exp(— iqr) , (8) 

after some calculation, we obtain the impact factors relevant for the total cross section: 



She N, 



2N,N^- 



- j dz£r\^{z, r) |2 (f((1 - z)qr) + F{zqr) - -^F{qr)^ . (9) 



The light-cone wave function for the g — ?■ QQ transition can be obtained from the well- 
known case for the photon as p, |lO[ : 






(^2 + (1 - zf^mlK^irnqr) + mlKl{mqi 



(10) 

where Kqi are generalized Bessel functions, and in the spirit of collinear factorization, 
we took the gluon to be on-shell. 



B. Dipole-dipole cross section 

It is now convenient to introduce the total cross section for two color dipoles of sizes 
r, s, in the two-gluon exchange Born-approximation |8|: 



aMr,s) =^^^^ I r„, ' F{qr)F{-qs 



4^^^«|(x(0) - xir) - xis) + xir - s)) , (11) 



where 



/72 
r 2 I ^212 exp(igr) = 4- if^cr) Kiificr) ■ (12) 

[Q + f^Gi f'G 

The Born level dipole-dipole cross section now reads 

(JDDir,s) = —^ ^ 1 - ncrKidxcr) - ij,gsKi{hgs) + ficlr - s\Ki{nG\r - s\ 

(13) 



N^ ^j^l 



Notice, that it is finite for /ic — )■ 0. 

The total cross section for the parton-level process gg — )■ QQQQ can now be written 
in terms of the dipole-dipole cross section and the light-cone wave-functions for the 
g — )■ QQ transitions as 



(^tot= I dzd^rdud^s\'^{z,r)\^\^{u,s)YT.{z,r]u,s) . (14) 

where 

■|az)z)((l-2;)r, (1-m)s) + aDD{{^ - z)r,us) - —aDD{{^ - z)r,s) 
+aDDizr, (1 - u)s) + aDoizr, us) - —aDoizr, s) 

c 

-j^[(yDD{r, (1 - u)s) + (TDD{r, us) - j^(yDD{r, s)j I . (15) 



C. Momentum distributions 

The mixed representation given above could in principle be used to obtain distribu- 
tions in longitudinal momenta, simply by stripping off the dz and/or du integrations. 
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FIG. 2: Feynman diagrams for the g — )• QQ impact factor. 

For the more interesting transverse momentum distributions it is better to start all over 
in momentum space, where the impact factor reads 



^\q,gn 



(fk r 1 

dz-:—^l [^{z, fc - (1 - z)q) - ^{z, k + zq)] —== 5ab\[QQ]i; z, k) 



(2vr)^ 



V2K 



+ [^{Z, fc - (1 - Z)q) - ^{Z, k + zq)] - dabc\[QQ]8, z, k) 

+ [^{z,k - {1 - z)q) - ^{z,k + zq) -2^{z,k)]^zfatc\[QQrs;z,k)y (16) 



After squaring, we obtain naturally the same structure as in |6[: 

I'^iz, k, q) = ^{^^^ [|vI/(^, k) - ^{z, k + zq)\' + \^iz, k) - ^{z, k - (l - z)q)\' 

-|^(^, A; - (1 - z)q) - ^(;z, k + zq)\^ + |^(2, fc - (1 - z)q) - ^(2;, k + zq)\'^\ . 

Sbc 



2N, 



I{z,k,q) 



(17) 



The explicit form for the squares of the light-cone wave functions can been found in [6| 
and reads: 



^(z,p) -^(z,p + k) 



as { [z^ + (1 - zf] 



P 



p + n 



p2 + mi, (p + k) 



m; 



+m; 



'^^p^ + rn^Q {p + K)'^ + mQ 



:i8) 



With all of this given, we can put together the differential cross section 
da 



N^-1 47i^al rr , ^rr , ^ , ^''k ^ dH d^ 
I[z,k,q)I[u,l,—q)dz .„ du- 



N'c [q' + ^^l] 



212 



'(27r)2 (27r)2 (27r)2 • 



(19) 



A brief comment on the kinematics of quarks is in order. The four-momenta of quarks are 
fully specified by the variables 2, A;, g, m, I. Incoming gluons carry longitudinal momentum 
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fractions a;i,X2 of the momenta of the incoming protons. The latter are 

Pa = P^n^ = y|n+, P, = p-n- = y^n" , (20) 

so that 2Pa ■ Pb = s. We parametrize four-momenta in hght-cone coordinates 

P={P+,P-,P), 2p+p^ - p^ = m^ . (21) 

Then, the four momentum, say of the quark/antiquark belonging to parent gluon 1, read 

We still need to give the relation between the transverse momenta Pq q of the (anti- 
)quark and the momenta fc, q used above. These relations read: 

PQ = k + zq, pq = -k + {l-z)q, (23) 

or, alternatively 

q = PQ + PQ, fc = (1 - z)pq - zpQ . (24) 

Notice, that q is the total transverse momentum of the QQ-pair, while k is the light-cone 
relative transverse momentum which was conjugate to the dipole size. Now having the 
four momenta, we can calculate all sorts of kinematical variables, e.g. the rapidity will 
be given by 

y« = iiog(|). (25) 

The kinematics of the "lower" QQ-pair is treated analogously. Obviously the trans- 
verse momentum of the second pair is just —q, because incoming gluons are collinear. 
When constructing four-momenta, we only need to be careful that the large component 
is along the lightcone-minus direction 

And here 

Pq = 1- uq, Pq = -1-{1- u)q . (27) 

D. On normalization and phase-space 

It is important to remember, that the expressions we derived are valid in a high energy 
limit, in which the invariant masses of QQ-pairs are much smaller than the center-of-mass 
energy squared of the gluon-gluon collision: M^2) ^34 ^ ^- 



In the practical calculation of the hadron-level cross section, we integrate over all 
momentum fractions carried by initial state gluons, and hence over the cms-energy of 
the gluon-gluon subprocess. In practice, we want that our cross section behaves smoothly 
also at low energies. 

Firstly notice, that our states are normalized in such a way, that the formulas are 
simple in the high energy limit. In particular, the constrained two-body phase space is 
just dz(Pk/{2'Kf. 

Effectively, the four-body phase space in the high-energy limit is just 

cPq , (fk , dH 

^ dz-——du-—-r. (28) 



(27r)2 (27r)2 (27r)2 ' 

Obviously it doesn't vanish if s approaches the threshold. To improve upon this, let us 
first rescale the amplitude, and introduce the Feynman-amplitude Aip 



(29) 



Mf = 2P+ V-2(l - z){2n) 2P_ Vm(1 - u){27v) ■ A 
= v^i(r^^^)2(2^y^I(r^^I)2(2^ -s-A, 

which is normalized such, that the fully differential cross section takes the form: 

da = — \MF\'^d^{s;pi,p2,P3,P4)- (30) 

2s 

Here, the constrained n-body phase space, with P^^ = s is 

d^s-p„ ...,pn) = (2vr)^<5W(P,„ - ^p,) . J] ^ 6{p^ - ml) . (31) 

We now introduce invariant masses M12, M34, and write the four-body phase-space as 

dM'^ dM"^ 

d<l>{s;p^,p2,P3,P4) = d^s;Pi2,P34) ^d^M^^;p,,P2) ^d^M^^;ps,P4) . (32) 

Now, going over to light-cone coordinates, we have 

•^'^'ta.P,-P.)n||;|. (33) 

Here we integrated out the Pj"-components from the on-shell conditions. Let us now 
write Pi = zP^, then one of the overall delta-functions gives us, that p^ = (1 — z)P^. 
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Furthermore, we can use the on-shell condition P^ 
then we obtain finally 



M^ to write 2P7oP^o = M^ + P 



12- 



27r 



d^{M'-pi,P2) 



dz d^p^d'^p2 ^(2)/ , o N 



(iz 



d'^k 



z{l - z) 2(27r)3 



(34) 



Which agrees, modulo the factors now absorbed in the Feynman amplitude with the 
result obtained in the high-energy limit. The major simplification in the high-energy 
limit occurs in the first factor of the four-body phase space (!32l) . Namely we can write 
the phase space for two clusters M^2: ^Ia ^ ^- which are seperated by a large rapidity 
distance as 



dHs; P„ P,) = {2nr 6^'\P - P, - P,) H ' ^^;^, ' Si2PtPr 



M^) . (35) 



For brevity, we wrote Pi = Pu, P2 = -P34 Now we can neglect the plus- component of P2 
and the minus-component of Pi in the overall four-momentum conservation, so that 



5(4) (p _p^_ p^) = 5(p+ _ p+)S{p- - P2)5^^\Pi 



(36) 



Hence all the integrations except for the transverse momentum ones can be done imme- 
diately. From the integrals over ±-components, we only get a factor of 



1 



2P+2P2- 



1 

2i 



(37) 



so that, finally: 



d^{s-Pi,P2) = - 



1 d^Pid^P2 



(27r) 



5'^'\Pi + P2)^i, 



d^q 



2s (27r)2 



(3J 



If we approach the threshold, the exact phase-space goes to zero and our approximation 
is very bad. Still, we know, that the integrated two-cluster phase-space will be just 



d^s-Pi,P2 



v/[s-(Mi2 + M34)2][S-(Mi2 



Mu? 



Svrs 



(39) 



In order to improve our calculation, we therefore introduce the correction factor 

/corr ~ 



1- 



(M12 + M;,Y 



1- 



(M12 - M34) 



2n 



(40) 



which within the approximations of the high-energy limit of course is exactly unity. A 
deviation from unity, or as a matter of fact any energy- dependence of the subprocess 
cross section thus indicates for us how far we are from the high-energy domain. 



E. pp — )• {QQ){QQ) inclusive cross section 

The cross section for proton-proton (proton-antiproton) can be calculated as usually 
in the parton model as 



^PP^iQQKQQ) 



(W) = / dXidX2 g{Xi,fJ.l) ^(x2,/i|) (^gg^{QQ)iQQ){s^^^) , (41) 



where s = XiX2W^. The last factor is the elementary cross section discussed in the 
previous sections. In calculating ( 14T1) we take into account kinematical constraints and 
the threshold correction factor ( HOl) . The elementary cross section is calculated first 
on a subsystem energy grid and a simple interpolation is done then when using it in 
formula (HT]) . We shall use different parton (gluon) distributions from the literature. 
The factorization scale of the gluon distribution in principle depends on the kinematics 
of the final state quarks. We use /i|, = 4mQ when calculating the integral fHT]) . The 
parton formula (HT]) is very useful to make differential distribution in invariant mass of 
the {QQ){QQ) system. In principle it will be desirable to obtain results for open-charm 
mesons, but the inclusion of hadronization goes beyond the scope of the present paper 
where we wish to present only a first estimation of the cross section for SPS production 
of {cc){cc). 



III. RESULTS 

A. gg -^ icc){cc) 

To make final calculations we have to fix a^ in formulae (fTSjl and ( IT9|) . We shall 
use leading-order running strong coupling constant as{kf + ml) in impact factors (see 
Eq. (TT8|) ) and o.s{[J^,) with ji^ = m.ax{kf^ + mg, q^) for gluon exchange (see Eq. (TT9|) ). 

In Figj3]we show total cross section for the gg — )■ {cc){cc) as a function of gluon-gluon 
energy. We show cross section with extra cut M12 + M34 < W and with extra correction 
factor (see Eq. (HCT]) ). The latter cross section will be used then to calculate corresponding 
cross section for proton-proton collisions. In the following calculations we have fixed the 
regularization parameter ^g = 0.5 GeV. There is only a marginal dependence on the 
value of the nonperturbative parameter. 

Before we go to real observables let us show an auxiliary distributions in qt and fcj 
(see Fig.gD). 

The transverse momentum distributions of quark (antiquark) is shown in Figj5] for 
a few selected subsystem energies. The higher subsystem energy the bigger transverse 
momenta are available kinematically. 

Rapidity distributions of quark (antiquark) for different energies are shown in FigEi 
We show distributions for quarks emitted from the upper line (solid line) and from 
the lower line (dashed line). The higher the energy the better separation of the two 
contributions can be seen. 

The rapidity separation can be better seen in the distributions in rapidity distance 
between (anti) quark- (anti) quark. The distance between quark-antiquark from the same 
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FIG. 3: Energy dependence of the total cross section for the gg — )■ {c.c)(cc) process. The 
description of the soHd hne is given in the main text. 
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FIG. 4: Auxihary distributions for the gg — t- {cc){cc) process for W = 20, 50, 100, 200 GeV 
(succesively growing cross section). 

pair is very small compared to the distance bewtween (anti) quark- (anti) quark from dif- 
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FIG. 5: The quark (antiquark) transverse momentum distribution da/dpt in the gg — >• {cc){cc) 
for W = 20, 50, 100, 200 GeV. 

ferent pairs. In FigJT] we show an example for subsystem energy W = 200 GeV. The 
distance in rapidity between quarks and antiquarks emitted from different pairs reminds 
a bit situation in double-parton scattering. There the distance between quarks and 
antiquarks emitted from two different hard processes can also be large |3|. 

Finally we close presentation of our results for gg — )■ (cc) (cc) by showing distributions 
in quark- antiquark invariant masses. Here there are two distinct classes of subsystems. 



We shall introduce the notation: 



da 
dMi2 



da 
dM-n 



da' 

dMi. 



(category I) 



for quark-aniquark emitted in the same pair and 

(category II) 



da 
dMi3 



da 
dMi4 



da 
dM23 



da 

dM2A 



da' 



dMj. 



for quark-antiquark emitted from different pairs. One can see that the average invariant 
mass of quark-antiquark from the same pair is smaller than the average invariant mass 
from the different pairs. At large invariant masses the emission from different pairs 
dominates over the emission from the same pair. The situation reminds that for double- 
parton scattering [3|]. 



B. 



VV 



ccccX 



Let us come now to proton-proton scattering. In FigjO] we show distribution in ,^i = 
logiQ(xi) (^2 = logio(a^2)) for I^ = 7 TeV (LHC). We see that typical xi and X2 are not 
too small, of the order of 10~^ - 10~^. This is the region where the gluon distributions 
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FIG. 6: Rapidity distribution of quark (antiquark) in the gg — >• {cc){cc) for W = 200 GeV. The 
solid line is for emissions from the upper pair, the dashed line is for emissions from the lower 
pair. 



are relatively well known. There is no strong dependence of the cross section on the 
choice of gluon distribution function (GDF). 

The distribution in invariant mass of the {cc) {cc) system (equal to subsystem energy) 
is shown in FigJTOl We show distribution for CTEQ6 GDFs [11]. For comparison we 
show distribution obtained for double-parton scattering (see |3|). While the double 
parton scattering contribution dominates in the region of small invariant masses, the 
single parton scattering contribution takes over above M^^ > 500 GeV. This is the 
region where large rapidity gaps between quarks and antiquarks occur. 

The energy dependence of the pp — )■ ccccX inclusive cross section is shown in FigJTTl 
One can observe that the inclusive cross section for the cccc final state is much smaller 
than that for the cc final state but grows somewhat faster at low energies. At higher 
energies the ratio is almost constant of the order of 1%. This is in contrast to double 
parton scattering contribution for the cccc production which grows much faster than the 
cross section for single cc production [3] . 



IV. CONCLUSIONS 

We have presented for the first time formulae for the production of two pairs of QQ in 
single-parton scattering. The elementary gg — ;■ {QQ){QQ) cross section was given in two 
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FIG. 7: Distribution in rapidity difference between quark-antiquark from the same pair (dashed 
line) and quark-quark, quark-antiquark, antiquark-quark, antiquark-antiquark from different 
pairs in the gg — t- {cc){cc) for the subsystem energy W = 200 GeV. 



different representations: so-called mixed one (longitudinal momentum fraction, impact 
parameter) called also dipole representation, and momentum space one within a high 
(subsystem) energy approximation. While the dipole representation is easy to include 
energy dependence of the dipole-dipole interaction the momentum representation seems 
better suited to include threshold effects. We have discussed how to correct the high- 
energy formulae close to threshold where the phase space is rather limited by energy- 
momentum conservation. 

We have presented energy dependence as well as different differential distributions for 
the elementary cross section for gg — )• {cc){cc). We have shown that the elementary cross 
section varies quickly from the kinematical threshold [W = Arric) up to W^ = 100 GeV 
where almost a plateau can be observed. The invariant mass distributions -^j^- = -^^^ 
are much steeper than that for :j^^ = -j^^ = -:;^— = -^^. where 1 and 2 are from the 
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first cc pair and 3 and 4 are from the other cc pair. 

Our high-energy approach does not include processes when a second pair is produced 
via gluon emitted from quark/ ant iquark of a first pair and its "subsequent" splitting 
into cc. This processes may be important for small rapidity distance between quarks 
and antiquarks. A good example can be LHCb kinematics when both c quarks (or both 
c antiquarks) are emitted within two units of rapidity. This mechanism should not be, 
however, important for the case of large rapidity interval between cc or cc discussed 
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FIG. 8: Quark-antiquark invariant mass distributions in the gg — )■ {cc){cc). The category I 
(same pair emission) distributions are shown by the dashed hue and category II (different pair 
emission) distributions by the solid Une. 

recently in the context of double parton scattering [3\. 

The elementary cross section has been convoluted next with gluon distributions in 
the proton. We have presented inclusive cross section for pp — >■ {cc){cc)X as a function 
of incident energy as well as invariant mass distribution of the cccc system. The results 
have been compared with corresponding contribution of double-parton scattering. We 
have found that the single-parton scattering contribution is significantly smaller than 
that for the double-parton scattering. We conclude therefore that a measurement of two 
pairs of cc at LHC would be very useful in testing models of double parton scattering. 

An evalution of distributions for charmed mesons would be very useful in planning 
and interpreting current measurements at the LHC. The LHCb collaboration started 
already such an analysis [12]. 
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FIG. 9: Distribution of the proton-proton cross section in ^i (^2) for the production of {cc){cc) 
at W = 7 TeV. The sohd hue is for fx'jp = m^ and the dashed line is for /i|, = Am^. In this 
calculation a simple leading order gluon distribution from [11] was used. 




s = 7 TeV 



1000 



M„ 



(GeV) 



FIG. 10: Distribution of the pp — )• ccccX cross section in M/^c for the production of {cc){cc) 
at for VF = 7 TeV. The solid line is for /x|n = ml and the dashed line is for /x|, = 4m^. In this 
calculation a simple leading order gluon distribution from [ll[] was used. 
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FIG. 11: Energy dependence of the cross section for the pp — )• ccccX reaction compared with 
that for the pp — )■ ccX reaction. 
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